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We present and study a novel electroweak-scale leptogenesis mechanism occurring in U{1)f flavor 
symmetric two Higgs doublet models with Higgs-dependent Yukawa couplings. In this scenario CP- 
violation originates entirely from the Higgs sector, and large CP asymmetries in TeV scale heavy 
neutrino decays can be obtained without resonance enhancement and any fine tuning of model 
parameters. Distinctive predictions for fermion Yukawa couphngs together with CP-violating Higgs 
sector make tests of the electroweak-scale leptogenesis mechanism realistic at LHC experiments. 
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INTRODUCTION 

The origin of quark and lepton Yukawa couplings re- 
mains the major unsolved question in particle physics 
One commonly accepted approach, the Froggatt-Niclsen 
mechanism explains the observed fermion mass hi- 
erarchies with powers of flavour symmetry U{1) p break- 
ing vacuum expectation value (VEV), ((S'f)/Af)", where 
Af is the scale of new physics, and n is determined by 
fermion quantum numbers under U{\)f- An interesting 
recent proposal [1, 01 suggests that all Yukawa couplings 
are generated by the SM Higgs boson H via the operators 
Cn{H^ H / M^)"' . This scenario predicts the flavor physics 
scale M - 0(1) TeV. However, as WH is a U{1)f sin- 
glet, no Froggatt-Nielsen type models can be built with 
one H, and the Yukawa hierarchies are based on the ad- 
ditional assumption that some coefflcients c„ vanish or 
on some unknown non-Abelian symmetry [3]. 

Leptogenesis Q is the most promising mechanism to 
explain the observed Q baryon asymmetry of the Uni- 
verse. In the standard scenario [5|] the CP asymmetry 
e^V; in heavy Majorana neutrino Ni decays is induced by 
complex neutrino Yukawa couplings y'l^j. In the case of 
0(1) TeV scale heavy neutrinos the seesaw mechanism 
imphes y"^ < 10"'', and large enough CP asymmetry 
cannot be generated unless masses of the heavy neutrinos 
satisfy the resonance condition , which requires unnat- 
ural fine tuning between Mpf. at low scale. 

In this Letter we propose a model for Higgs-dependent 
Yukawa couplings which successfully explains the masses 
and mixings of quarks, charged leptons and neutrinos. 
We extend the SM Higgs sector with another doublet. Be- 
cause HuHd is nontrivial under U(1)fj Froggatt-Nielsen 
type model building defines Yukawa coupling hierarchies 
from fermion fiavor quantum numbers [j]. In two Higgs 
doublet models the CP-violation may come from the 
Higgs sector [§] . We show that in this scenario new lepto- 
genesis mechanism occurs below the scale of electroweak 
symmetry breaking (EWSB), T < Tc, in which the nec- 
essary CP violation comes entirely from the Higgs poten- 
tial. The CP asymmetry is not suppressed by small neu- 
trino Yukawa couplings and successful leptogenesis oc- 
curs naturally for non-degenerate heavy neutrino masses 



independently of the initial conditions on Ni abundances. 
Taking into account flavor effects in all heavy neu- 
trino decays we derive and solve Boltzmann equations 
for the lepton asymmetry as well as for sphaleron tran- 
sitions [llj. Although sphalerons decay soon below the 
EWSB scale, < [l3, , we show that the lepton 
asymmetry is safely converted to the baryon asymmetry. 
This scenario implies a phenomcnological upper bound on 
Mn- from successful leptogenesis. Because the Higgs sec- 
tor CP-violation is testable at the LHC Q , our proposal 
opens completely new, realistic, perspectives for test- 
ing electroweak-scale leptogenesis mechanisms at collider 
experiments. 

THE MODEL 

We consider global U{1)f symmetric two Higgs dou- 
blet model with flavor charges given by Q{Hu) = 
hu, Q{Hd) = hd, QiU) = -4K - 3hd and Q(iVfl,) = 0. 
We assume hu + hd ^ so that the term HuHd is 
not U{1)f invariant. The U{1)f invariant Yukawa La- 
grangian for heavy Majorana neutrinos Nfh is given by 



HuHd 
AP 



1 



NrM'^N'r + h.c, (1) 



and similarly for quarks and charged leptons. In Eq. ([T|) 
we have neglected the term Nj^L^Hl (HuHd/M^)'''''^^ 
because its contribution to the induced Yukawa couplings 
is additionally suppressed. Eq. H]) gives rise to {2n'^ + 1) 
Higgs boson interactions with N and L. Below the EWSB 
scale the Higgs bosons acquire VEVs (Hu) = Vu = w sin/3 
and {Hd) = = ucos/3, where v = 174GeV. The small 
flavor symmetry violating parameter in this scenario is 
e = VuVd/M^ ~ 10~^, and for tan/3 = 1 the cut-off 
scale of this model is M = 1.23 TeV. Expanding the 
neutral components of Higgs fields {v + Hq)^"' "'"-'^/M^" 
in Eq. ([1]), we get the fermion mass terms proportional 
toe" f and multi-Higgs interaction terms e" Hq{Hq/vY , 
k — 0, .., 2ny , which have the common suppression factor 
e" . The existence of quadratic interaction e" H^Nv/v 
below EWSB plays a crucial role in our scenario. 

Under the charge assignment presented above, n^^ are 



universally given by 



3. This gives the correct 



2 



pattern for neutrino Dirac Yukawa couplings and, to- 
gether with the Majorana mass terms, imply the seesaw 



N 



induced light neutrino mass scale rrii, 
1.5 X 10~^ eV for — ITeV. As is common to Abelian 
flavor models, the observed neutrino masses and mix- 
ing are obtained by appropriate adjustment of the 0(1) 
coefficients j/^-*^. In the charged lepton sector the flavor 
charges are given by nfj = 3, n^j 



2 and n|^- 



imply the measured masses via rur ~ 



1, which 



and me ~ 3e vcp . Similarly, the correct quark masses and 
mixing are obtained with flavor charges nf^ 



^2 



f3 = l,n«. = (2,2,l),n-.=n-. = (1,1,0). 
The Higgs potential of this model is given by 

y _ ™2 ,2 



\Hu 



\Hdf + M\Huf + >^2\Hdf 



+ A3|-ff„|'|i?d|' 



A4 \HuHd\ 



+ 



m^HuHd + A5 {HuHdf + Ae \Hu\^ HuHd 



+ Ay \Hd\ HuHd 



(2) 



(3) 



where Eq. ^\sU{1)f symmetric and Eq. ([3]) is explicitly 
U{\)f breaking part of V. While m\ ^ and Ai^2,3,4 are 
real, m? and As^e,? are in general complex and give rise to 
CP-violation. Minimization of the potential V requires 

h = m^~ v'^spcp (A4 + 2A5 - Ag tan /3 - Ay cot (3) (4) 

to be real and the vacuum stability condition implies 



Ai > 0, A2 > 0, A3 + VA1A2 > 0, 



A3 + A4 + 2VX^- 2IA5I > (for Ae = Ay = 0). (5) 

The mass matrix {1/2)HqMqHq for neutral Higgs bosons 
Hq = (Rei/°, ReiJ)^, ImiJ", Im_ff°)^ has one zero eigen- 
value which corresponds to the Goldstone boson eaten by 
Z. Mq is diagonalized by the orthogonal transformation 
O^MqO — dia.g{mf^,mf^^,m\^,0), and the mass eigen- 
states ha, a — 1, ..,4, are obtained from the weak eigen- 
states Ho via h = O^Hq. Since CP is violated in the 
Higgs potential by complex parameters, real and imagi- 
nary components of the neutral Higgs bosons are mixed 
with each other. If all couplings were real, off-diagonal 
blocks of Mq vanish, and hi^2 and /13 would correspond 
to the CP even and odd Higgs bosons, respectively. 

Inserting Higgs mass eigenstates into Eq. ([1]) implies 
^B. ^ interactions up to dimension-five operators as 



A%K + -Bt^Kht, 

J V 



(6) 



where the effective Yukawa couplings A and B are 



r>ab 



V2 '^'^ " 

<■ (02a + iOia)] , 



{Ola+iOsa) 



(7) 



(-1)^ 



X [nj;- - 3) (O16 + zOsb) 

+ <j «i - 1) (02fc + iOib) C/3] (a • 



C0 



b). 



(8) 



The couplings A and B are complex because O as well as 
are complex. However, the CP asymmetry in heavy 
neutrino decays depends on \y''\^ and, therefore, only the 
phases in O contribute to the CP asymmetry. 
The Ne. — gl interactions are given by 



Ce = (-l)<:'e"-yr,7V.PLe,0+ + h.c. 



(9) 



where </)"'" is the charged Higgs boson with mass 2h/ sin 2/3. 
The three-point vertex for neutral Higgs bosons is 



V3 = w Cabchahhc, 



(10) 



where the CP-violating couplings Cabc are lengthy ex- 
pressions of the parameters in V. 

LEPTOGENESIS 

In our scenario the heavy neutrino decays N — > LH 
are induced below EWSB. The usual leptogenesis CP 
asymmetry induced by Ni loops Q is negligible for non- 
degenerate 0(1) TeV neutrinos. However, in neutral de- 
cay modes Ni — > X]a=i ^j^^a new CP asymmetry 

is generated by interactions ([6]) and (|lGp via the diagrams 
in FIG. [H In Eq. pT|l 
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where Ja,bc is the loop function defined by 
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Ja,bc = / dx In [x(x - 1) + (1 - x)rb + xr^ 



(12) 
(13) 

(14) 



where , 



/m\ . The function (fH)) has absorptive 



imaginary part when the mass of final state ha is larger 
than the sum of intermediate Higgs masses, nih^ +mh^ < 
m^^. These states induce the CP asymmetry PT|) . 

To exemplify the magnitude of generated CP asym- 
metry we assume non-degenerate heavy neutrino masses 
{MN)ij = MA,diag(0.5, 1.25, 1.5), Mn = 1 TeV, and 



Ai 

™2 



0.2, A2 0.5, ImAg ^ 0, 

(300GeV)^ + 2iv^spcplm.\z, others = 0, (15) 



and tan /? = 1 . The imaginary part of rrP' is constrained 
by the minimization condition Eq. For this choice 
the CP violating coupling ImAs is the only free parame- 
ter. The neutral Higgs boson masses are approximately 
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FIG. 1; Diagrams generating dominant contribution to CP 
asymmetry in heavy Majorana neutrino decay below EWSB. 
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FIG. 2: Examples of CP asymmetries as functions of ImAs. 
Decays N3 — > '^^i'2ha give the largest contribution to e. 
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FIG. 3: Evolution of |77s-L/£max| and rjNi with z = Mni/T 
for ImAs ~ 10~^ in the case of thermal (solid curves) and zero 
(dashed curves) Ni initial abundances. 



drj. 



3 



dz 



eq 



VAel e] 



(18) 



where t^aj and the asymmetry in left-handed neutrinos 
(charged leptons) ?7A!/(e)j are related to each other by the 
"A-matrix" defined as »7Ai/(e)j — ^'ji^'^VAk with 



hi - 140GeV, /i2 ^ 480GeV, hs ~ 440GeV and de- 
pend weakly on the magnitude of ImA5. Therefore, the 
Ni decay modes to /i2,3 (with hi in the loop) have the 
absorptive part and non-vanishing CP asymmetry. 

Three components of nine CP asymmetries are plot- 
ted in FIG. [Has functions of ImAs. The shaded region is 
prohibited by the vacuum stabilization condition Eq. ([5]). 
For the chosen parameters the decays N3 J^a ^^ha 
give the dominant contribution to e. Therefore we denote 
e| = Emax- As sceu in FIG. [2J large CP asymmetries of 
orded 10~^ can be generated without resonant condition 
for heavy neutrinos. The dependence of e on A5 is linear. 

In addition, the heavy neutrinos Ni also decay via 
Eq. (HI) with the width Ff,^^ = F(A'^^ -> el^cj)^), 



Wikl 



(16) 



No CP asymmetry is generated in these decays. The total 
decay rate of Ni is thus given by Toi = ^'hi^^'hi^ where 

'- Di — l^k=l ^ Di "^^^^ ^ Di — 2^k=l ^ Di- 

In order to calculate the generated baryon asymmetry 
we derive Boltzmann equations for the evolution of heavy 
neutrino and Aj = B/3 — Lj asymmetry abundances rj^i 
and 77Aj , respectively. Here rjx = nx/n^, where nx is 
the number density of X and the photon number density 
is — 2T^ (,{?>)/ Ti"^ . The Boltzmann equations read 



dz 



njH{z = 1) 
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iDi 



(17) 
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(19) 



(20) 



The thermally averaged decay rates jDi^'oi'' ^^'^ I'oi''^ 



are obtained from F/^i, F^'^'^'' and F^'^'^''"' in Eqs.((T3 



by 7Di 



and 



^j = l IDi 

modified Bessel functions and z. 



E 



nZ{Ki{z,)/K2{z.))rjjf', Inf = 
and 7£)j = 7]^^ -|- -yf^^ . Ki and K2 are the 
Mmi/T. Hereafter we 
define z = zi as usual. 

FIG. [3] shows the evolution of Ifys-i/smaxl and rusn 
with z for the case of thermal (solid curves) and vanish- 
ing (dashed curves) Ni initial abundances. In both cases 
realistic B — L asymmetry tib-l ~ Ej "HAj is generated 
below the critical temperature Tc, z^ = Afjvi/Tc — 3.04. 

Fast sphaleron processes convert the generated lep- 
ton asymmetry to baryon asymmetry [ll|. How- 
ever, sphalerons decay quickly at Td < Tc IJ, 
Il3| . The sphaleron rate Fa(s+l) at temperatures 
MwiT) < T < Mw{T)/aw is given hyTAm+L) ~ 
Mw{Mw/awT)^{Mw/Tfexp[-Esp/T] 0, [l| where 
aw is the SU{2)l fine structure constant, is W- 

boson mass and the sphaleron energy is Esp ~ Mw/aw 
Just below Tc the sphaleron interactions are faster than 
the expansion rate of the Universe, ^a{b+l)/H{z = 
1) ^ 1. In this region, rjB-L is converted to baryon 
and lepton asymmetry 77s, l by sphaleron effect with the 
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FIG. 4; Lepton and baryon asymmetries for ImAs = 10~^ 
(blue and green curves) and ImAs — 10"'' (purple and red 
curves) for thermal {th) and vanishing (0) initial abundances 
of Ni ai z = Zc- The horizontal band is the experimental value 
r?|j"P = (6.05 - 6.37) x 10-^° [|. 



successful leptogenesis implies an upper bound Mpf < 4.5 
TeV on the heavy neutrino masses. The distinct predic- 
tion 0, 01 of factor of three (five) enhancement of t, b (/i) 
Higgs-dependent Yukawa coupUngs compared to the SM 
prediction is, in our model, affected only by a small cor- 
rection related to the generated baryon asymmetry (by 
As effects in our numerical example). While discover- 
ing TeV scale Ni with small Yukawa couplings at LHC 
is more than a challenging task [l^, this prediction, to- 
gether with CP-violating Higgs sector testable at LHC 
[9[, makes our scenario of electroweak-scale leptogenesis 
realistically testable at LHC. 
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temperature-depedent rate given by [1^, |l£ 
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30T2 + 21w(T)2 
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(21) 
(22) 



The sphaleron rate rA(B+L) decreases below Tc 
by the Boltzmann factor cxp[—Esp/T] and reaches 
^A{B+L){zd)/H{z = 1) — 1 a,t Zd — 6.04. Since the 
sphaleron processes are effectively switched off at z > z^;, 
the baryon asymmetry is unaffected below z^j. 

Numerical results for the evolution of baryon and lep- 
ton asymmetries are shown in FIG. [3] for ImAs ~ 10^^ 
(blue and green curves) and ImAs = lO"'^ (red and pur- 
ple curves). Indeed, the sphalerons decouple at Zd while 
the lepton asymmetry still evolves at z > z^. For appro- 
priate As the observed baryon asymmetry can be easily 
generated both for thermal and vanishing initial Ni abun- 
dances denoted by (th) and (0), respectively. 

Because the leptogenesis window between Zc and Zd 
is fixed, successful leptogenesis implies an upper bound 
Mjv < 4.5 TeV on the heavy neutrino mass scale. For 
larger masses Ni decay above the EWSB scale without 
being able to generate the observed baryon asymmetry. 

CONCLUSIONS 

We have presented a novel mechanism for electroweak- 
scale leptogenesis in the two Higgs doublet scenario of 
Higgs-dependent Yukawa couphngs with global U{1)f 
flavor symmetry. In our mechanism the necessary CP- 
violation comes entirely from the Higgs potential, and 
large CP asymmetries can be naturally obtained with- 
out resonance enhancement from heavy neutrinos nor any 
other fine tuning of model parameters. There is a small 
window of temperatures below Tc where sphalerons are 
active, yet the observed baryon asymmetry can be eas- 
ily obtained before sphalerons decay. The requirement of 
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